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SUMMARY 


General parameter-plane equations are derived to generate stability bound- 
aries for a class of systems characterized by a feedback loop that contains a 
complementary filter and a model for either the low- or high-frequency portion 
of the plant. This combination allows those frequencies of the part of the 
plant that is modeled to be fed back for control while suppressing other 
frequencies. 

For all specific examples considered, the stability regions obtained using 
the complementary filter and frequency model were larger (and in some cases, 
considerably larger) than those obtained using a low-pass filter in the feed- 
back of the system output. Furthermore, higher gain control was possible. 


INTRODUCTION 

A filter can be thought of as a circuit, a differential equation, or a 
transfer function which passes certain frequencies in an input while attenuat- 
ing or rejecting others. The specific filter discussed in this paper is 
referred to as a complementary filter. It consists of two subfilters - a 
high-pass filter and a corresponding low-pass filter - the outputs of which 
are summed. 

A primary application of the complementary filter is as follows. Suppose 
two different sensors are used to measure a state of a system. Moreover, let 
the first instrument be good at measuring the high-frequency content of the 
signal but not the low-frequency content. In other words, the first instrument 
records the signal plus a low-frequency noise. The second instrument gives an 
accurate indication of the signal at the low frequencies but degrades at the 
higher frequencies (high-frequency noise) . One approach is to pass the signal 
from the first instrument through a high-pass filter and the signal from the 
second instrument through a low-pass filter and then sum the results to get a 
representation of the whole signal. The only constraint placed on this system 
by the complementary filter is that the transfer functions of the high-pass and 
the low-pass filters must sum to unity. Under ideal conditions of no noise, 
a signal which passes through a complementary filter will be reproduced exactly. 
No distortion or phase shift will occur in the signal. However, with noise, 
this will not be the case. With a noisy input, parameters in the complementary 
filter are adjusted so that the total impurity from the low-frequency noise in 
the first instrument and the high-frequency noise in the second instrument is 
minimized. 

The complementary filter is a simple concept, but at the same time, it is 
a very powerful and useful concept. The complementary filter in earlier work 
(refs. 1 and 2) is referred to as a "distortionless" filter. According to the 
editorial comment in reference 1 , this type of filter enjoyed much success in 
a variety of applications around 1956. It has applications in aircraft flight 



measurements and inertial navigation (refs. 3, 4, and 5). A compensation- 
filtering scheme was used in reference 3 to significantly increase the band- 
width of an early Skylab control-moment-gyro system configuration by combin- 
ing noisy position- sensor information and "clean" rate-command information. 

For practical reasons, the rate command was used instead of the actual system 
rate. In this sense, the compensation filter is a complementary filter. Brown 
(ref. 4) discussed the relationship of the complementary filter to the Kalman 
filter. Higgins (ref. 6) continued this comparison and showed that the comple- 
mentary filter is a special case of a steady-state Kalman filter. Schmidt 
(ref. 7) compared the complementary and Kalman filters for the STOLAND system. 

It is important to note at this point that application of the complementary 
filter does not require any knowledge of the more complex Kalman filter. 

Accurate mathematical models of the high-frequency modes of a system are 
not always available. Consequently, in a control design based on models of the 
low-frequency modes, it may be undesirable to feed back the very high frequen- 
cies, which were not taken into account in the control design. Also, in flexi- 
ble booster control, it is important not to set up any resonance with the struc- 
tural mode of the system. At the same time, however, it is desired to control 
the rigid-body modes. A novel and ingenious use of the complementary filter 
was made by Tutt and Waymeyer (ref. 8) in flexible booster control: "This 

approach does not adapt to body bending, but instead is contrived to ignore 
it . . ."by using the basic concept of the complementary filter and a rigid- 
body model in the feedback loop. This same basic* idea was recently applied in 
the high-vibration environment of a helicopter to eliminate low-frequency rotor- 
induced vibrations (refs. 9 and 10). This latter work used a high-frequency 
plant model. 

The objective of this paper is to apply the parameter -plane method (ref. 11) 
to derive stability boundary equations for a particular class of systems 
(defined by a general block diagram) , which use the Tutt and Waymeyer (ref. 8) 
feedback approach. The parameter-plane method is used to display the stability 
region in the two-dimensional space of a parameter of the complementary filter 
and a control gain. Specific examples are presented, and a standard low-pass 
filter feedback is used for comparison. 


SYMBOLS 

A (I , J) matrix element in the Ith row and Jth column 

A function of S defined in deriving explicit expressions for K 

and T 

Ar,Ai real and imaginary parts of A, respectively 

B function of S defined in deriving explicit expressions for K 

and T 

Br,Bi real and imaginary parts of B, respectively 
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C R' C I 

D 

D R' D I 
D (I) 

D(S) 

D(S) 

F (S) 
G(I) 

H (S) 

h 

i 

I,J 

K 

Ki,K 2 

l 

M 

N(I) 

N(S) 

N(S) 

S 

1st 

T 

T 1 rT 2 


function of S defined in deriving explicit expressions for K 
and T 

real and imaginary parts of C, respectively 

function of S defined in deriving explicit expressions for K 
and T 

real and imaginary parts of D, respectively 
denominator of the Ith transfer function G(I) 
denominator of system transfer function H(S) 

characteristic equation, or denominator of H(S) relatively prime 
to N (S) 

Laplace transform of input signal to complementary filter 
transfer function of the Ith system 
system transfer function 
high-frequency noise 

= R 

integers 

system control gain constant 

particular values of K 

low-frequency noise 

function defined in equation (37) 

numerator of the Ith transfer function G(I) 

numerator of system transfer function H(S) 

numerator of H(S) relatively prime to D(S) 

Laplace complex variable 
absolute value of S 

parameter in low-pass and high-pass filters 
particular values of T 
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t time 

U(S) Laplace transform of system input 

X L' X H outputs of low-frequency and high-frequency plant models, respectively 
(Xg) Xjj after passing through low-pass filter 

X(S) Laplace transform of system output 

Y(S) Laplace transform of output signal from complementary filter 

a, functions defined in equations (33), (34), and (35), respectively 

X, p quantities which take on the values of either unity or negative unity 

o real part of S 

co imaginary part of S 


ANALYSIS 

The complementary filter is briefly discussed, then the system block diagram 
for the present analysis is presented, followed by a formulation of the charac- 
teristic equation of the system. Using this characteristic equation for the 
block diagram, the parameter-plane method (ref. 11) is used to develop general 
equations for generating stability boundaries in the plane defined by a filter 
parameter and the system control gain. Stable regions are identified from these 
stability boundaries. 


Complementary Filter 

Figure 1 shows the complementary filter operating on an uncorrupted signal 
F (S) and having an output Y(S). Since H(S) and 1 - H(S) occur in parallel, 
the resultant transfer function is given by their sum, which is unity - hence, 
the name complementary filter. The point to be made from figure 1 is this: 

If a signal F(S) passes through a high-pass filter with transfer function 
H (S) and a low-pass filter with transfer function 1 - H(S) , which are comple- 
mentary, and if the signals from the two filters are summed, then the output 
signal Y(S) will be the same as the input signal F(S). There is no distor- 
tion or phase shift of the input. 

Figure 2 shows a complementary filter with noisy inputs. Suppose Z is 
a low-frequency noise and h is a high-frequency noise. The output is 
described by 
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Y(S) = H(S)[F(S) + A] + h - H(S)Nf(S) + h] 
= F(S) + {H(S)£ + [l - H(S)]h} 


( 1 ) 



Notice in equation (1) that Y(S) = F(S) 
in the braces. 

Figure 3 shows the specific form of 
in this paper. Equation (1) becomes 


Y(S) = F (S) + — V + (_L 

\TS + 1/ \TS 


except for the residual error term 
the complementary filter of interest 



In signal processing, the problem is to choose T to minimize the residual 
error. The following section uses the same filter components shown in figure 3 
but with different inputs to accomplish a desired objective. The influence of 
T on system stability is emphasized. 


Frequency Feedback Model With Complementary Filter 

Figure 4 shows the system block diagram examined in this paper. The system 
has an input U(S) and an output X(S). The overall plant is composed of low- 
frequency and high-frequency plants preceded by a servo. The complementary 
filter consists of the specific high-pass and low-pass filters shown in fig- 
ure 3. The objective is to examine stability regions in the plane of the filter 
parameter T and the system gain K. 

The operation depicted in figure 4 is explained as follows: Let X L and 

X H be the outputs of the low-frequency plant and high-frequency plant, respec- 
tively. If the low-frequency plant model and servo model are exact models, then 
the output of the low-frequency plant model will also be X L . Notice that the 
input to the high-pass filter component of the complementary filter is X L , while 
X L + Xjj is the input to the low-pass filter component. Since X L passes 
through both the high-pass filter and low-pass filter components of the comple- 
mentary filter, it will be unchanged. Meanwhile, X H passes through the low- 
pass filter. Hence, as shown, the resultant feedback is X L + ( x h)lP' that is, 
the output of the low-frequency plant plus the output from the high-frequency 
plant which has been passed through a low-pass filter. The term ( x h)lp 
the low-frequency part of the high-frequency plant. Adjusting T adjusts the 
low-frequency part of X H that is fed back for control. This paper examines 
the effect of K and T on system stability. 


Characteristic Equation 

The notation in figure 4 is convenient for deriving general expressions 
involving the system transfer functions. The N(I)'s and D(I)*s are the 
numerators and denominators, respectively, of the block transfer functions. 

The system transfer function between X(S) and U(S) in figure 4 is 
denoted by 
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1 


X(S) 

U(S) 


= H (S) 


(3) 


with 


N(S) N(S) 

H(S) = = (4) 

D(S) D (S) 


where N(S) and D(S) signify polynomials in S which may have some common 
factors that can be canceled. The numerator N(S) of the transfer function 
and the denominator D(S) have no common factors and are mathematically said 
to be relatively prime polynomials in S. 

The polynomial D(S) is referred to as the characteristic polynomial, and 
D(S) =0 is called the characteristic equation. Roots of D(S) are called 
characteristic roots. It is well known that a linear, time- invariant system 
is stable (X ( t) •*0 as t °°) if and only if all roots of the characteristic 
equation have negative real parts. 

The system transfer function H(S), expressed in terms of the individual 
G (I) block transfer functions in figure 4, is easily shown to be 


K[G(1) + G ( 2 ) ]G (6) 

H(s) = (5) 

1 + k[G(4) G (5) G (7) + G (1 ) G(3) G (6) + G(2) G(3) G (6) ] 


Replacing the G(I)'s in equation (5) with their equivalent ratios shown in fig- 
ure 4 allows equation (5) to be written as equation (4) with 


N(S) = K[N(1) D (2) + D (1 ) N(2)]N(6) D(1) D(2) D(3) D(4) D<5) D(6) D(7) (6) 


D (S) = D 0 ) D(2) D(6){D(1) D(2) D(3) D(4) D(5) D(6) D(7) 
+ k[D(1) D (2) D (3) N{4) N(5) D(6) N(7) 


+ 

N(l) 

D(2) 

N{3) 

D(4) 

D (5) 

N(6) 

D(7) 


DO) 

N(2) 

N(3) 

D(4) 

D (5) 

N(6) 

D(7) ]} 

Notice in figure 4 

that 
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N(4) = TS 


(8) 


D (4) 

= D (3) 

(9) 

N(3) 

= 1 

(10) 

D (3) = TS + 1 

perfect modeling, G(5) = G(l), and G(7) = G(6); therefore, 

(11) 

N(5) 

= N (1 ) 

(12) 

D(5) 

= D(l) 

(13) 

N(7) 

= N (6) 

(14) 

D(7) 

- D(6) 

(15) 


For convenience, equation (6) is expressed as 

N(S) = D (1 ) D(2) D (4) D (5) D(6) D(7) N(S) (16) 

where 

N(S) = k[n(1) D (2) + D (1 ) N(2)]D(3) N(6) (17) 

Equations (9), (13), and (15) allow equation (7) to be factored as 

D(S) = D (1 ) D (2) D (4) D (5) D(6) D(7) D(S) (18) 

where 

D(S) = D (1 ) D (2) D (3) D (6) + k[d(2) N(4) N(5) N(7) 

+ N (1 ) D (2) N (3) N (6) +D(1) N(2) N(3) N (6) ] (19) 
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Ill 


The ratio of N(S) to D(S) is the system transfer function. Recall that 
D(S) =0 is the characteristic equation. By equations (8), (10), (11), (12), 
and (14), the characteristic equation can be expressed as 

D (1 ) D (2) (TS + 1) D (6) + K[D(2) (TS) N(l) N(6) + N(l) D(2) N(6) 

+ D (1 ) N(2) N(6) ] = 0 (20) 


A more convenient and compact form of the characteristic equation for appli- 
cation of the parameter-plane method is 


TA + KTB + KD + C = 0 (21 ) 

where 

A = S[(D(1) D (2) D (6) ] (22) 

B = S [D (2) N (1 ) N (6) ] (23) 

C = D(1 ) D (2) D (6) (24) 

D = N (1 ) D (2) N (6) + D (1 ) N(2) N(6) (25) 


The values of S which satisfy equation (20) or equation (21) are the char- 
acteristic roots. Clearly, these roots will vary as K and T take on differ- 
ent values. The parameter-plane method is used to examine how the stability con- 
dition (stable or unstable) changes with different combinations of K and T. 


Stability Boundaries 

In this section, equations are derived to generate stability boundaries 
in the plane of K and T for the system shown in figure 4. For comparative 
purposes, similar equations are developed for a low-pass filter feedback system. 
This latter system is the same as the system in figure 4, except that the path 
through the servo model, frequency model, and high-pass filter is removed. 

Complementary filter system .- The numerators and denominators of the block 
transfer functions in figure 4 are functions of S. Consequently, A, B, C, 
and D in equations (22) to (25) are functions of S. Let S = iu> and write 
A, B, C, and D as complex numbers: 


8 


I 


A = Ar + 

iAi 

(26) 

B = B R + 

iBi 

(27) 

+ 

II 

a 

iCi 

(28) 

D = Dr + 

iDi 

(29) 


Substituting equations (26) to (29) into equation (21) and setting the 
resulting real and imaginary parts equal to zero yields the two simultaneous 
equations for K and T: 

TA r + K (TB r + D r ) + C R = 0 (30) 

TAj + K (TBj + D x ) + Ci = 0 (31) 

Subtract equation (31) multiplied by (TB R + D R ) from equation (30) multiplied 
by (TBj + D-j ) to get the quadratic equation in Ts 

OtT 2 + BT + Y = 0 (32) 

where 


a = BjA r - BrAj 

(33) 

3 = DjAr + BjCr - DrAj - BrCj 

(34) 

Y = DjCr - DrCj 

(35) 

The roots of equation (32) are 


-3 ± 


T - 

(36) 

2a 
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M = 3 2 - 4ay 


(37) 


The two roots in equation (36) are denoted as T-| for the positive radical and 
as T 2 for the negative radical. The gain K, obtained by adding equations (30) 
and (31), is 


T(Ar + Aj) + Cr + Cj 

K (38) 

T(Br + Bj) + Dr + Dj 


In equation (38), K = K-| when T = T-| , and K = K 2 when T = T 2 . 

At this point, it may be worthwhile to state explicitly the meaning of equa- 
tions (36) and (38) . For particular constant values of K and T, the charac- 
teristic equation will have a set of corresponding roots. As K and T are 
varied, these roots will move, generating root-locus curves. Suppose it is 
desired to find values for K and T which result in a root-locus curve inter- 
secting the imaginary axis at some specified point S = iu)*, where 10 * is a 
particular value of w. This is accomplished by setting S = iio* in the char- 
acteristic equation, equating the real and imaginary parts to zero, and solving 
the two resulting simultaneous equations for K and T. 

Equations (36) and (38) produce those combination values of K and T 
which will result in a root-locus curve intersecting the imaginary axis at 
S = iw. For a given value of to, the value of T is calculated from equa- 
tion (36) , which is an implicit function of a). Then, K is calculated from 
equation (38) using this value of T and the same value of w. By letting 
to vary, a plot of K versus T can be generated. The resulting curves are 
stability boundaries that correspond to the imaginary axis in the S-plane. 

Hence, in the same manner that the imaginary axis divides the S-plane into a 
stable region (left half-plane) and an unstable region (right half-plane), the 
stability boundaries partition the plane of K and T into stable and unstable 
regions. 

Since only real values of T are of practical interest, calculations are 
only continued for M = 0 in equation (37). Also, since the characteristic 
roots occur in complex-conjugate pairs, it is sufficient to plot T versus K 
as u) increases from zero. Using negative values of to would only duplicate 
the results obtained for positive values of (0. 

Low -pass filter system .- The low-pass filter system is the same as the sys- 
tem shown in figure 4, except that the path passing through the servo model, 
frequency model, and high-pass filter is omitted. It can be shown that the 
associated characteristic equation is 
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D(S) = D (1 ) D (2) D (3) D(6) + k[n(1) D(2) N(3) N(6) 
+ D (1 ) N (2) N(3) N(6) ] = 0 


(39) 


Equations (10) and (11) permit equation (39) to be written as 


D (1 ) D (2) (TS + 1) D (6) + K[N(1) D(2) N(6) + D(l) N(2) N(6)] = 0 (40) 


or, in the more convenient form. 


A (TS +1) + KB = 0 


(41) 


where 


A = D (1 ) D (2) D (6) (42) 

B = [ N (1 ) D (2) + N (2) D (1 )] N(6) (43) 

Setting S = iu) and using the complex representations for A and B 
(eqs. (26) and (27)) modifies equation (41) to 


(Ar + iAj) Tia) + (Ar + iAj) + K (Br + iBj) — 0 


(44) 


Setting the real and imaginary parts in equation (44) equal to zero results in 
the two simultaneous equations: 


-AjUa) + Ar + KBr = 0 


ArTdo + Aj + KB j = 0 


The simultaneous solutions of equations (45) and (46) are 


A r Bi - B r Aj 

T = 

0)(AiBi + B r A r ) 


- (a r 2 + Ai 2) 

K = — ^ 

ArBr + AjBj 


(45) 

(46) 


(47) 


(48) 


In equations (47) and (48), T and K are single-valued functions of u. 
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Zero characteristic roots .- It is only natural that in solving for T to 
obtain equations (36) and (47) , division by zero is encountered. This follows 
because T always appears with S in the transfer functions for the low-pass 
and high-pass filters in figure 4. With S = iw, one can think about solving 
for ujt and then dividing by w to get an equation for T. The u is shown 
explicitly in the denominator of equation (47) for emphasis. Hence, in generat- 
ing the stability boundaries, special attention must be given to the question: 
What combination values of T and K result in the characteristic equation 
having zero roots (S = 0) ? The computer program which formed the system charac- 
teristic equation from the component transfer functions and then solved for its 
roots used the following convenient notation. 

Let 


N(l) 

= A(1 ,1) 

+ A(1 ,2) S 

+ A(1 ,3)S 2 + . . . 

(49) 

D(l) 

= A(2,l ) 

+ A(2,2)S 

+ A(2,3)S 2 + . . . 

(50) 


and so forth. Then, form the matrix where the first row contains the coeffi- 
cients of the polynomial N(1), the second row the coefficients of the poly- 


(1 ) , etc. Symbolically, 



N(l) 

A (1 ,1) 

A (1 ,2) 

A(1 

/ 3) 

D(l) 

A(2, 1 ) 

A(2,2) 

A(2 

,3) 

N(2) 

A(3,l) 

• • • 



D (2) 

A(4,l ) 

• • • 



N (3) 

A(5,l) 




D (3) 

A(6,l) 




N(4) 

A(7, 1 ) 




D (4) 

A(8,l) 




N(5) 

A(9,l) 




D (5) 

A(1 0 , 1 ) 




N(6) 

A (11 ,1 ) 




D(6) 

A(1 2,1 ) 




N(7) 

A(1 3,1 ) 

• • • 



D (7) 

A (1 4,1 ) 

• « • 




( 51 ) 
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The constant term for N(l) is A(l r l); for D(1), is A(2,l), etc. Using 
this notation and setting S = 0 in equation (20) gives 


-A(2,l) A(4,1) A(1 2,1 ) 

K = (52) 

A(l,1) A(4,l) A(1 1 ,1 ) + A(2 , 1 ) A(3,l) A(1 1 ,1 ) 


For later applications, N(6) = 1 or A(ll,l) = 1 , so that equation (52) 
reduces to 


-A (2,1 ) A(4 , 1 ) A(1 2,1 ) 

K = 

A ( 1 ,1) A(4,l ) + A(2,l ) A(3,l) 


(53) 


When K has the value computed using equation (53) , the characteristic 
equation has a root on the imaginary axis at the origin, regardless of the value 
of T. This constant value of K, for arbitrary value of T, is plotted also 
in the parameter plane as a stability boundary. In the K,T plot, this special 
boundary appears as a vertical line which intersects the K-axis at the value 
assigned by equation (53). 

Infinite characteristic roots .- To obtain the stability boundaries associ- 
ated with infinite characteristic roots, the combination values of T and K 
which satisfy the characteristic equation when |s| = | iw| 05 are computed. 
This is accomplished by dividing the characteristic equation by its maximum 
power in S. Afterwards, the only term not having some power of S in the 
denominator is the original coefficient of the maximum power of S. Hence, as 
| S | 00 , this coefficient approaches zero in order to satisfy the remaining 

characteristic equation. 

For all subsequent applications, the maximum power of S in the charac- 
teristic equation will occur in the D (1 ) D(2) (TS) D(6) term in equations (20) 
and (40) . After dividing by the maximum power of S and afterwards letting 
| S | °°, there remains a constant term times T equated to zero. After divid- 

ing by this constant term, T = 0. The result is that T = 0 and K arbitrary 
is a stability boundary, that is, the K-axis. 


Stability Regions 

As the coordinating parameter u > 0 increases in the set of equations 
developed for K and T, a plot of T versus K can be drawn as, for example, 
in figure 5. This partitions the T,K plane into different regions which are 
either stable or unstable. The stability condition of a given region is then 
evaluated by choosing a point inside the region and computing the characteris- 
tic roots for this point. The number of roots with positive real parts 
(unstable roots) in each region is circled in figure 5. A circled zero indi- 
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cates a stable region. The number of unstable roots gained or lost in moving 
from one region to another is evident. In subsequent applications, only the 
stable regions are emphasized. 


Calculation Procedure 

Figure 6 illustrates the general procedure used to construct the stability 
regions for the complementary- filter system. First, a set of transfer functions 
for the blocks in figure 4 is input to a computer program. Next, the quantities 
w, M, T-j , K] , T 2 , and K 2 are obtained from a computer printout generated 
by incrementing w. From this printout, values of u> for which M = 0 are 
selected for further consideration. The program has options for plotting 
T-j versus K-j or T 2 versus K 2 . To generate the stability boundaries 
(T versus K curves) , oj is incremented in the range where M = 0, and values 
of T (eq. (36)) are plotted against corresponding values of K (eq. (38)). 
These boundaries can be overlooked if this incrementation interval for W is 
too large. Efficient generation of the stability boundaries requires close 
monitoring of the computer printout and plotter by the investigator. 

The same procedure is followed for the low-pass filter system. The program 
has an option for selecting either the complementary or low-pass filter system. 
For the low-pass filter, equations (47) and (48) are pertinent. 


Degree of Stability 

In this paper, the substitution S = ito is used to generate stability 
boundaries. For values of K and T falling within the stable regions, con- 
fined by these boundaries, there are no characteristic roots located on or to 
the right of the imaginary axis in the S-plane. In a similar manner, stability 
regions could be generated such that for combination values of K and T fall- 
ing within these regions, there are no roots any closer to the imaginary axis 
than some specified number a. This is accomplished by setting S = cr + iu) and 
following the same procedure described in this paper. For further reading on 
relative stability, refer to references 11 and 12. 


High-Frequency Plant Model 

A low-frequency plant model is indicated in figure 4 because the objective 
was ultimately, through use of the complementary filter, to filter out the high 
frequencies from the feedback signal. 

In reference 9, it was desired to filter out the low frequencies; hence, 
a high-frequency plant model was used. To use the equations of this paper for 
a high-frequency plant model, simply call G(l) the high-frequency plant, 

G(2) the low-frequency plant, and G(5) the high-frequency plant model, and 
interchange G(3) and G(4). 
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Approximation Models 


All numerical results presented in this paper are based on perfect models 
of the servo and low-frequency plant of the real system. To examine the effects 
of modeling errors, equations (12) to (15) should not be used. This translates 
into the fact that equation (21) still holds, except that 


A = SfD(l) D (2) D (5) D (6) D(7)] (54) 

B = S[D(1) D (2) N(5) D (6) N(7)] (55) 

C = D (1 ) D (2) D(5) D (6) D(7) (56) 

D = N (1 ) D (2) D (5) N (6) D(7) + D(1) N(2) D(5) N(6) D(7) (57) 


Equations (54) to (57) allow variations in the models G(7) and G(5) from 
those of G (6) and 6(1), respectively. Everything else remains the same. 
There are no changes in the equations corresponding to the low- pass filter 
system. 

The flight investigation of reference 9, which used a high-frequency plant 
model, reported that the "... overall system performance was insensitive to 
mismatch between the dynamics of the plant and the plant model used in the com- 
plementary filter computation.” 


APPLICATIONS 

The overall plant is composed of two subplants, a high-frequency plant and 
a low-frequency plant. Arbitrary plant transfer functions are assigned for numer- 
ical computation. For the high-frequency plant. 


N (2) 0.3S 2 + 0.45S + O.S 

D'(2) s 2 + 0.06S + 9 


(58) 


and, for the low-frequency plant: 


N (1 ) 2S + 1 

0(1) S 2 + 1 . 4XS + p 


(59) 


where X and y are either 1 or - 1 • 
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Four special cases are enumerated based on the denominator term in the low- 
frequency plant: 


Case 

1 : 

S 2 

+ 1 .4S + 1 

(stable plant: 

X = 1 , 

y - i) 

Case 

2: 

S 2 

- 1 . 4S + 1 

(unstable plant: 

X = -1 , 

M - 1) 

Case 

3: 

S 2 

+ 1 .4S - 1 

(unstable plant: 

X = 1, 

y = -i) 

Case 

4: 

s 2 

- 1 . 4S - 1 

(unstable plant: 

X = -1 , 

y = -i ) 


Stability regions resulting from feedback implementing a low-pass filter 
are compared graphically with those from feedback incorporating the complemen- 
tary filter in figure 4. Plots are presented for two conditions: (1 ) without 

servo (figs. 7 to 10) and (2) with servo (figs. 11 to 14). 


Without Servo 

Case 1 .- Figure 7 shows the block diagram and stable regions resulting from 
use of the low-pass filter and the complementary filter. Comparing figures 7(b) 
and 7(c) shows the improvement brought about by using the complementary filter. 
Especially note that larger K values are permitted. In effect, larger K 
values mean tighter control because K enhances the system error so that it 
can be further reduced. 

The stability boundaries in figure 7(b) and 7(c) intersect the K-axis at 
the same points. This will always be the case, because for T = 0 in figure 4, 
the complementary filter reduced to the low-pass filter. 

Case 2 .- The low-frequency plant has negative damping and is unstable. 

The block diagram and stability regions are shown in figures 8. Note the negli- 
gible stable region for the low-pass filter in figure 8(b) and the vast improve- 
ment in figure 8(c) for the complementary filter. 

Case 3 .- The results in figure 9 show a larger stable region for the com- 
plementary filter. 

Case 4 .- Figure 10 displays a tremendous improvement in the size of the 
stability boundary when the complementary filter is used over that when the low- 
pass filter is used. 


With Servo 

The block diagrams in figures 11, 12, 13, and 14 are the same as those in 
figures 7, 8, 9, and 10, respectively, except for the servo and servo model. 
The same basic conclusions persist, namely, (1) the system with the complemen- 
tary filter has a larger stable region in the K,T plane than the system with 
only a low-pass filter and (2) for a given value of T, a stable system can be 
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maintained while using larger error gains K for tighter control. Inciden- 
tally, figure 11(b) corresponds to the stable region shown in figure 5 for 
T = 0. 


CONCLUDING REMARKS 

A complementary filter is commonly used in signal processing to estimate 
a state or signal when more than one measurement of the state is available. In 
previous work, the complementary filter has been used in conjunction with a 
model of the rigid- body modes for flexible booster control. This approach 
allows the structural modes to be filtered from the feedback signal, while the 
low-frequency rigid-body modes are not. 

In this paper, general parameter-plane equations are derived to generate 
stability boundaries for a class of systems characterized by a feedback loop 
that contains a complementary filter and a model for either the low- or high- 
frequency portion of the plant. This combination allows those frequencies of 
the part of the plant that is modeled to be fed back for control while suppress- 
ing other frequencies. The plant and model rational transfer functions can be 
specified arbitrarily. 

Numerical examples are presented which compare the stability regions 
obtained using the complementary-filter-model system with those regions incurred 
in using just a standard low- pass filter feedback. For the selected examples, 
it is shown that the complementary-filter-model system produces a larger stabil- 
ity region and allows larger error gains (tighter control or regulation) to be 
used. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 9, 1980 
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Figure 1,- Complementary filter with noise-free input signal. 
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Figure 2.- Complementary filter with noisy inputs. 







































Specify system transfer functions (figure 4) 


Obtain vaiues for: oo, M, Tj, Ky T^, 
(Eqs. (36), (37), (38)) 


Plot stability boundaries (T versus K) 
(Eqs. (36), (38)) 


Identify stable regions by computing characteristic roots 

(Eq. (20)) 


Figure 6.- Calculation procedure. 
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S 4 + 40S + 400 


LOW-FREQ. PLANT 


2S + 1 


GH-FREQ. PLANT 


.3S 2 + .45S + .9 



S 2 + .06S + 9 


LOW-FREQ. 

SERVO MODEL 

PLANT MODEL 

400 

2S + 1 

S 2 + 40 S +400 

S 2 - 1.4S - 1 


HIGH- PASS FILTER 


TS 


TS + 1 


LOW- PASS FILTER 


TS + 1 

(a) Block diagram. 

Figure 1 4.- Block diagram and stable regions for case 4 with servo. 




















1. Report No. 

NASA TP-1744 

4 . Title and Subtitle 


2. Government Accession No. 


STABILITY BOUNDARIES FOR SYSTEMS WITH FREQUENCY- 
MODEL FEEDBACK AND COMPLEMENTARY FILTER 

7. Author(s) 


L. Keith Barker 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

15. Supplementary Notes 


3. Recipient's Catalog No. 
5. Report Date 

November 1980 


6. Performing Organization Code 
8. Performing Organization Report No. 

L-l 3787 

10. Work Unit No. 

506-54-13-02 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Paper 

14. Sponsoring Agency Code 


16. Abstract 

General parameter -plane equations are derived to generate stability boundaries for 
a class of systems characterized by a feedback loop that contains a complementary 
filter and a model for either the low- or high-frequency portion of the plant. 

This combination allows those frequencies of the part of the plant that is modeled 
to be fed back for control while suppressing other frequencies. For all specific 
examples considered, the stability regions obtained using the complementary filter 
and frequency model were larger (and in some cases, considerably larger) than 
those obtained using a low-pass filter in the feedback of the system output. Fur- 
thermore, higher gain control was possible. 


17. Key Words (Suggested by Author(s)) 


18. Distribution Statement 


Complementary filter 
Parameter plane 
Stability boundaries 
Frequency model feedback 

High-gain control _ 

19. Security Oassif. (of this report) T 20. Security Classif. (of this page) 


Unclassified - Unlimited 


Subject Category „66 

21. No. of Pages 1 22. Price 


Unclassified 


I Unclassif ied 1 40 I 

For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley , 1980 



National Aeronautics and 
Space Administration 


THIRD-CLASS BULK RATE 


Washington, D.C. 
20546 


Official Business 

Penalty for Private Use, $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 



10 1 1U,G, 102780 S00903DS 

DEPT OF THE AIR FORCE 
AF WEAPONS LABORATORY 
ATTN: TECHNICAL LIBRARY (SOL) 
KIRTLAND AFB NH 87117 




POSTMASTER: ^ Undeliverable (Section 158 

Postal Manual) Do Not Return 


